Thermalization and temperature distribution in a driven ion chain 
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We study thermalization and non-equilibrium dynamics in a dissipative quantum many-body 
system — a chain of ions with two points of the chain driven by thermal bath under different 
temperature. Instead of a simple linear temperature gradient as one expects from the classical heat 
diffusion process, the temperature distribution in the ion chain shows surprisingly rich patterns, 
which depend on the ion coupling rate to the bath, the location of the driven ions, and the dissipation 
rates of the other ions in the chain. Through simulation of the temperature evolution, we show that 
these unusual temperature distribution patterns in the ion chain can be quantitatively tested in 
experiments within a realistic time scale. 

PACS numbers: 



I. INTRODUCTION 



Many-body non-equilibrium dynamics has raised sig- 
nificant interest in recent years, in particular with 
connection to the atomic experiments where far-from- 
equilibrium phenomena can be conveniently investigated 
due to the long relaxation time of these systems P, Q . 
For instance, the dynamics after a quantum quench has 
been studied in several model Hamiltonians [3 6] . The 
physics becomes even richer if we add engineered dis- 
sipation to the underlying system, which is possible to 
realize in atomic experiments. Several interesting effects 
have been analyzed recently from interplay between in- 
teraction and dissipation in cold atomic gas [7|, 1|| . 

Motivated by this line of research, in this paper we 
study thermalization and temperature distribution in a 
linear ion chain, with two points of the chain driven re- 
spectively by a heating and a cooling thermal bath. In 
a sense, the configuration here is an analog of a classi- 
cal example in thermodynamics — the heat propagation 
in a bar with its two ends fixed at different tempera- 
ture. In contrast to a simple linear temperature gradient 
as one sees in a classical bar, we find that the temper- 
ature distribution in a driven ion chain shows very rich 
behaviors in its steady states: first, the temperature dis- 
tribution depends critically on the ratio between the driv- 
ing speed and the interaction rate in the ion chain. In 
the weak driving region, the temperature distribution is 
non-monotonic across the ion chain and shows a strange 
mirror effect. In the strong driving region, instead of a 
linear temperature distribution, all the ions between the 
two driven ones are stabilized to a constant temperature 
which is in the middle of the two bath temperatures. Sec- 
ond, we find that the bath which drives the ion with a 
medium coupling rate plays a more effective role in deter- 
mining the temperature of the other ions in the system. 
When the bath drives the ion strongly, it has little influ- 
ence on temperature of the other ions in the chain, which 
is pretty counter-intuitive. To see the unusual phenom- 
ena predicted in this paper, we discuss the requirements 
for an experimental observation and calculate the time 
dynamics to achieve the steady-state temperature dis- 



tribution. The qualitative features of the steady-state 
temperature distribution patterns are pretty insensitive 
to the size of the ion crystal and they show up already in 
a pretty small ion chain with less than ten ions. The real- 
ization time to the steady state is realistic for observation 
compared with the current experimental time scale. 

The paper is arranged as follows: In Sec. II we provide 
the theoretical model and the calculation method. The 
main results of the temperature distribution patterns un- 
der different circumstances are presented and discussed 
in Sec. III. In Sec. IV, we discuss the time dynamics to 
achieve the steady-state temperature distribution. Sec. 
V summarizes our major findings. 



II. MODEL 

We consider a chain of ions along the axial z-direction 
coupled with the Coulomb interaction and driven indi- 
vidually by a thermal bath. The thermal bath can be 
provided, for instance, by cooling or heating laser beams 
shined on each ion. Both the driving rate and the effec- 
tive bath temperature can be tuned by controlling the in- 
tensity and the detuning of these laser beams. The oscil- 
lation of the ith ion around its equilibrium position is de- 
scribed by the coordinate and the momentum operators 
Xi,pi, with the dynamics determined by the Heisenberg- 
Langevin equation: 



Pi 

Pi = ~^2jAijXj 



HPi 



(1) 



where Ay denotes the coupling matrix between the ions, 
7i is the driving rate of the bath, and Q(t) denotes the 
corresponding random force from the thermal bath. We 
consider in this paper the ions' motion along the trans- 
verse x direction, so the coupling matrix Ay is given by 

A ii = "1 ~ Ej(^i) 1/ \ z j - z if and A ij = V \ z j - z i\ 3 
for i 7^ j, where u) x is the transverse trapping frequency 
and Zi denotes the ion equilibrium position along the ax- 
ial direction. We take the ion spaci ng dp as the length 
unit, e 2 /do as the energy unit, and \J e 2 /(rndg) as the fre- 
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quency unit so that every quantities in the matrix Ay be- 
come dimensionless. Note that with control of an anhar- 
monic trapping potential along the axial direction, one 
can make the ion spacing uniform in a scalable trap Q. 
For the conventional harmonic axial trap, the ion spacing 
is not uniform. In that case, the equilibrium positions Zi 
are determined numerically with the given trapping po- 
tential, and we take the unit do as the largest spacing in 
the ion chain. The temperature distribution pattern that 
will be shown in Sec. Ill is insensitive to the details of the 
axial trapping potential, and for most of the calculations 
in the following, we assume the ion spacing is uniform for 
simplicity (with exceptions in Sec. Ill for consideration 
of the time dynamics in a harmonic trap). For inde- 
pendent Markovian bath, the random force Q(t) can be 

- b\) with the bosonic field 



expressed as Q = —i\/u)i/2(bi 
operator bi(t) satisfying (&j(ii)&j(i2)) 



where is the average phonon number which charac- 
terizes the temperature of the bath, and u>i = y/ An is 
the local oscillation frequency of the ith ion by fixing all 
the other ions in their equilibrium positions. For typical 
linear ion traps with strong transverse confinement, we 
have Ui &s uj x . The correlation of Q(t) is then given by 

(Q(hMh)) = "i( T f + i/2)M(*i - *a)- 

The Langevin equation (JTJ can be solved exactly 
through diagonalization, with the solution formally ex- 
pressed as 



q(i) 



-fit 



q(o) 



dTe n(T-t) »7(7 



(2) 



where q = a; 2 , Pi,Pz, — ) T 



{0} 



{Xi} 



, and f2 = 



is a 27V x 27V 



-I 

[Ay] [ 7i <5y 

matrix, which can be diagonalized as [C/ — 1 S7[/] a ^ — 
From this formal solution, we obtain the vari- 
ance of the operators Xi and Pi as 
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where [A = 1,2, ...,7V correspond to the x-operators and 
fi = TV + 1, TV + 2, 2N correspond to the p-operators. 
The temperature of each ion is similarly characterized by 
the average phonon number of its thermal oscillation, so 
the temperature of the ith ion as a function of time t is 
represented by T,(t) = (t)) + (p? l). As 

long as 7j is nonzero for some ions, each eigenvalue X a 
has a positive real part due to coupling of all the ions, 
and the system approaches to a steady state as t — 5- oo. 
The steady state temperature for each ion is denoted by 
T? = Ti(t — ► oo). 
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Figure 1: (Color online) The steady-state (t — >• oo) tempera- 
ture distribution Tf (measured by the mean thermal phonon 

number) of a uniform ion chain with TV = 100, where the l s ^ 
itlu 



and the last (100 u 
(at temperate Tj 8 



ions are driven respectively by a cooling 
2) and a heating (at temperate T 2 S — 10) 
thermal bath under varying driving rate 7. The temperature 
of the driven ions (l s ^ and 100^) are indicated by the thick 
curves. The other ions are assumed to be isolated from the 
bath, an d the tran sverse trapping frequency uj x =10 in the 
unit of a/ e 2 / (mdg), where do is the ion spacing. 



III. STEADY-STATE TEMPERATURE 
DISTRIBUTION 



In analogy to the example of heat propagation in a 
conducting bar, we consider a long ion chain with the 
two edge ions driven by different thermal bath with tem- 
perature and T 2 S (with Tf < T 2 S ), respectively. 
For simplicity, we assume the corresponding driving rate 
7i = 72 = 7 and all the ions in the middle have no dis- 
sipation. After all the ions attain the steady state, the 
temperature distribution Tf across the chain is shown in 
Fig. Q] for N = 100 ions under different driving rates 7. 
Apparently, the temperature distribution does not fol- 
low a linear gradient. There are three regions for the 
distribution of T?, depending on the ratio between the 
driving rate 7 and the ion interaction rate. Note that in 
our unit the interaction energy between the neighboring 
ions is of the order of unity. If 7 3> 1, the dissipation 
is much faster than the energy propagation in the chain, 
and without surprise the two edge ions have temperature 
basically fixed by their corresponding bath temperature 
Tf and T 2 B . However, it is surprising that all the other 
ions in the middle approach almost the same temperature 
given by Tf ~ (T/ 3 + T 2 s )/2 for i = 2, 3, • • • , 99 in this 
case. The temperature does not fall down gradually from 
the hot end to the cold end as in the classical heat prop- 
agation problem, but takes a sharp jump right from the 
driven ion to the next one and then keeps constant over 
the whole chain. In the opposite limit of weak driving 
with 7 <C 1, the phonon propagation is faster than the 
bath driving, and all the ions approach the same temper- 
ature. It sounds that the temperature distribution in this 
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Figure 2: (Color online) The steady-state temperature distri- 
bution for N = 101 ions, where the cooling bath (at T\ = 2) 
remains at the l s ^ ion and the heating bath (at = 10) 

moves to (a) the middle (51 s ^) ion, and (b) the 30^ ion. The 
temperature of the driven ions are indicated by thick (red) 
curves along the axis of the driving rate 7. The other pa- 
rameters are the same as in Fig. 1. The temperature profiles 
are highlighted by different thick (green) curves along the ion 
index axis to represent 7 = 10~ 3 , 10 -2 , 0.1, 1, 10. For small 
7, note that the temperature profiles both in (a) and (b) are 
found to have the mirror effect, i.e., a reflection symmetry of 
the distribution with respect to the middle (51 s ^) ion. 



limit resembles a classical thermal equilibrium, however, 
this picture is not true. The example in Fig. Q]represents 
an exception instead of a rule where we put the cooling 
and heating ions exactly at the symmetric positions of 
an ion chain. As we will see in the following, when we 
shift the position of one of the ions to break the reflec- 
tion symmetry, the temperature distribution in the weak 
driving limit has strange features. It is even not mono- 
tonic across the ion chain, in sharp contrast with the 
distribution from the diffusion process. Between these 
two limiting regions there lies a transition region with 
0.01 < 7 < 1, where the temperature of the edge ions 
gradually approach the corresponding bath temperature. 
Note that the temperature for the ions next to the driven 
ones follow a non-monotonic curve as one changes 7, al- 
though such variation is pretty small. 

To break the position symmetry of the cooling and the 
heating ions, we fix the cooling ion still at the edge, but 
move the heating ion inside the chain. Fig. Uta) shows an 
example of the temperature distribution where the heat- 
ing ion is right at the center of the chain. When 7> 1, 
the heating ion separates two plateaus in the tempera- 



ture distribution: a lower temperature region between 
the heating and the cooling ions and a higher temper- 
ature region on the other side of the heating ion. As 7 
decreases, the two plateaus smoothly descend to the same 
temperature while the heating ion remains as a peak in 
the temperature distribution. When 7 < 1, the temper- 
ature distribution has a reflection symmetry with respect 
to the chain center. A dip in temperature appears on the 
free edge of the chain (the left side of Fig. [Ha)), mir- 
roring the cooling ion on the other side of the chain (the 
mirror effect). The temperature distribution is appar- 
ently non-monotonic in this case. The mirror effect exists 
for other positions of the heating ion. For instance, Fig. 
&b) shows the temperature distribution where the dis- 
tance between the heating and the cooling ions is about 
one third of the chain length. Both the cooling and the 
heating ions have their mirror images in the temperature 
distribution in the weak driving limit. 

In the above calculation, we assumed 71 = 72 for the 
cooling and the heating ions. When these two driven 
ions are put at the edge, the middle ions are at a con- 
stant temperature T m which is exactly the average of 
the two bath temperatures. If the driving rate 71 7^ 72, 
we may expect that the bath which drives the ion more 
strongly plays a larger role in determining the temperate 
T m of the middle segment ions. This expectation, how- 
ever, turns out to be not true. To show that, we map 
out the middle segment temperature T m in Fig. [3] as a 
function of 71 and 72. The figure shows that 7$ 0.1 
is the optimal driving rate for which the corresponding 
bath has the largest influence on the temperature T m . At 
this optimal value, the driving rate is comparable with 
the ion interaction rate in term of the order of magni- 
tude. The middle segment temperature is close to the 
temperature of the bath that drives the ion at the opti- 
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Figure 3: The middle-ion steady-state temperature T„ for the 
configuration shown in Fig. 1, as a function of two indepen- 
dent bath driving rates 71 and 72. The dark (light) regions 
correspond to low (high) temperatures, and the numbers there 
indicate different temperature contours. The parameters are 
the same as in Fig. [T] expect that 71 associated with the cool- 
ing bath(at Tj 8 = 2) and 72 associated with the heating bath 
(at T 2 S = 10) are varied independently. 
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Figure 4: (Color online) (a) The steady-state temperature 
distribution across the ion chain, when the other ions in the 
chain (except the two driven ones at the edge) are coupled 
to a background thermal bath at temperature T® g = 4 with 
a coupling rate 72 = ... = 799 = jbg- (b) The temperature 
profiles at different 7^ (cross-section view of (a)). We find 
that the temperature profile for the main part of the chain, 
excluding the two edge ions, approaches a linear distribution 
at an intermediately small 7^ ~ 10~ 2 7. 



mal rate. When the driving gets too strong or too weak, 
the bath plays little role in determining the temperature 
of the other ions in the chain. This result has important 
implication for the sympathetic cooling [Iol - [l^ |: instead 
of fast cooling of the ancilla ions, cooling at a moderate 
optimal rate is more efficient in reducing the temperature 
of the computational ions. 

So far we have neglected dissipation of the middle seg- 
ment ions. Now, apart from the cooling and the heating 
bath (with temperature T-f and T%, respectively) at- 
tached to the two edge ions, we assume all the middle 
ions are coupled to a background bath with temperature 
T bg ( T F < T bg < T ?) at a coupling rate 7fcff . The final 
temperature distribution of the ions is shown in Fig. 0] 
under different background coupling rates 7t> 9 . In this 
calculation, we fix 71 = 72 = 7 = 0.1. For tiny 7b g , 
the temperature distribution shows no noticeable differ- 
ence compared with the case of jt>g = 0. However, as 
jbg increases to a moderate value with Jbg/j ~ 0.01, 
the temperature distribution of the middle segment ions 
has a clear linear spatial gradient, resembling the linear 
temperature distribution of a classical bar in the heat dif- 
fusion problem. As "fbg further increases and gets close 
to 7, the temperature of the middle ions are pined to the 

bg 



temperature gradient disappears again in this limit. 



IV. TEMPERATURE EVOLUTION 

We note that the temperature distributions shown in 
this paper are insensitive to the size of the ion system. 
We have checked the temperature distribution of the ions 
with the size of the chain varying from ten to a few hun- 
dreds of ions and noticed no qualitative change of the dis- 
tribution pattern. To experimentally test these unusual 
temperature distribution patterns, one can start with a 
small system of a few ions that are within the reach of 
the current experimental technology. The final tempera- 




background bath temperature T b B as one expects and the 
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Figure 5: (Color online) The temperature relaxation dynam- 
ics for an ion chain with N — 20 in (a) a uniform anharmonic 
trap, (b) a harmonic non-uniform trap, and (c) a harmonic 
trap with a background coupling rate 7^/7 = 10 -3 and a 
background temperature Tj% = 4. The 1 st ion and the 20 th - 
ions are driven respectively by a cooling (at T± = 2) and 
a heating (at T 2 S = 10) bath, with the same driving rate 
71 = 72 = 7 = 0.1, as indicated by the thick (red) curves. 
The initial temperature is assumed to be Ti(t = 0) = 5 for 
all the ions. The transverse frequency u x = 10 in the unit of 
sj e 2 /(mdg), where do is the largest spacing in the ion chain 

(between the 1 st and the 2 nd , or the 19 th and the 20 th ions) 
for the harmonic cases. In terms of real numbers, with a typ- 
ical choice of do = 10/im for 20 ytterbium ions in a harmonic 
trap, the trapping frequencies for the transverse and the ax- 
ial traps are given respectively by uj x as 2ir x 1.4MHz and 
u) z ~ 27T x 76.5kHz, and the driving rate 7 « 14kHz. 



5 



ture (the mean phonon number of the ion motion) can be 
detected, for instance, by measuring the scattering side- 
bands of a laser beam through a ion. The asymmetry in 
the blue and the red sidebands and their ratio give di- 
rect inference of the thermal phonon number of the ion 
motion [iH-[l5l|. 

To probe the properties of the steady states, we need 
to know the time scale to approach these steady states. 
Fig. [5] shows the relaxation dynamics for 20 ions in either 
a harmonic or an anharmonic uniform trap . For a uni- 
form trap (in Fig. [5^), the calculation shows that there 
are two time scales in the thermalization. First, with a 
time scale t\ ~ 1/7, the temperature of the two edge 
ions quickly approaches the corresponding bath temper- 
ature. During this step, the temperature of the middle 
segment ions only changes slightly, and the change gets 
smaller as one moves away from the driven ions. After 
that, a longer time scale ti sets in, representing the inter- 
action driven thermalization process. All the ions grad- 
ually approach the steady-state temperature. Note that 
the temperature of the edge ions (as a well as the other 
ions that are close to the two driven ones) does not fol- 
low a monotonic evolution curve. Instead, it first comes 
pretty close to the corresponding bath temperature and 
then is dragged back toward its steady state value. The 
value of the second time scale increases with the sys- 
tem size (roughly linearly) and t% ~ 40/7 for 20 ions. For 
a ytterbium ion ( 171 Yb + ) chain with spacing 10 /im, £2 
is about 3 ms, which is a pretty reasonable time scale for 
experiments. For a harmonic trap (in Fig. EJd), the basic 
feature is similar except that the second time scale t% gets 
site dependent. The edge ions and their neighbors ap- 
proach the corresponding steady-state temperature with 
a time scale that is comparable with the case of a uniform 
chain. However, as one moves away from the edge, the 
thermalization time gets much longer with an exponen- 
tial increase. For the middle ion, the thermalization is 
not finished yet with t ~ 10 10 /7. With such an extremely 
long time scale, of course one can not neglect the small 



dissipation of the other ions to the background bath. If 
we take into account a small background dissipation, e.g., 
with a rate "fbg/j ~ 10 -3 as shown in Fig. [5}:, the long 
tale in the temperature evolution is completely gone for 
the harmonic trap. Now, the temperature of all the ions 
approach their steady values within a time scale that is 
comparable with the uniform case, although the steady- 
state temperature of the middle segment ions is dragged 
down a bit toward the background bath temperature. 



V. CONCLUSION 

In summary, we have shown that the steady state tem- 
perature distribution of a driven ion chain shows surpris- 
ingly rich patterns. Many features of these patterns are 
unexpected, and they are in sharp contrast with the sim- 
ple linear temperature gradient as one sees in the classical 
heat diffusion problem. Our calculation is based on exact 
numerical methods with no unreliable approximations, 
so we believe all the unusual temperature distribution 
patterns revealed by this calculation will show up in ex- 
periments, although we still lack intuitive explanation of 
many of these features. We also investigate the relaxation 
dynamics and the time scale to reach the steady state, 
and show that these patterns should be observable under 
a realistic time scale in a small system that is within the 
current experimental reach. 
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